Quantum Magnetic Impurities in Magnetically Ordered Systems 
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We discuss the problem of a spin 1/2 impurity immersed in a spin S magnetically ordered 
background. We show that the problem maps onto a generalization of the dissipative two level system 
(DTLS) with two independent heat baths, associated with the Goldstone modes of the magnet, 
that couple to different components of the impurity spin operator. Using analytical perturbative 
renormalization group (RG) methods and accurate numerical renormalization group (NRG) we show 
that contrary to other dissipative models there is quantum frustration of decoherence and quasi- 
scaling even in the strong coupling regime. We make predictions for the behavior of the impurity 
magnetic susceptibility that can be measured in nuclear magnetic resonance (NMR) experiments. 
Our results may also have relevance to quantum computation. 
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Quantum impurity problems are characterized by a 
single quantum mechanical degree of freedom coupled to 
a reservoir. These are the simplest problems in physics 
that exhibit non-trivial many-body effects. Among them, 
the DTLS [1] plays a central role because it is related to 
a variety of different physical processes such as the dissi- 
pative tunneling of a particle in a double well or the cou- 
pling of an Ising spin to a gapless fermionic environment, 
i.e, the anisotropic Kondo problem. The Kondo prob- 
lem is one of the best understood impurity problems and 
has been studied by a large variety of methods. Its ther- 
modynamic properties can be studied exactly via Bethe 
ansatz [2] and RG [3] , its dynamic properties at low ener- 
gies can be calculated via conformal field theory [4], and 
many of its correlation functions can be obtained exactly 
[5]. The problem of quantum impurities immersed in 
magnetic media close to a quantum phase transition has 
also attracted a lot of attention recently due to its pos- 
sible relevance to cuprates, heavy-fermions and organic 
materials [6]. While most of the recent works focus on 
the paramagnetic phase we concentrate on magnetically 
ordered phases [7,8]. 

In this paper we study a quantum impurity problem 
of a different nature, namely, the problem of a single 
spin 1 /2 coupled to a c?-dimensional magnetically ordered 
system with spin S. A possible application of our re- 
sults can be found, for instance, in KMni_a;Cu2:F3 with 
X < 1 (Cu has spin 1/2 and Mn has S" = 5/2). KMnFa 
is a three dimensional (d=3) cubic quantum Heisenberg 
antiferromagnet (QHAF) with exchange coupling J be- 
tween the Mn spins. This material orders in a Neel state 
and has well defined gapless magnon modes [9] . The Mn 
spins interact with the Cu spin via an exchange J' . In 
the past this problem was studied by a series of differ- 
ent theoretical techniques that usually assume the impu- 
rity spin to be aligned with the surrounding Neel state 
and/or J' fa J [9,10]. As we are going to show, while 



this kind of approach is warranted at low temperatures, 
it fails to describe the full quantum behavior of the im- 
purity at energy scales intermediate between J' and J 
(J' ^ J) where the impurity spin fluctuates strongly 
and gets mixed with the quantum fluctuations of the 
magnetic environment. The problem at hand is simi- 
lar to the Kondo effect of magnetic impurities in metals 
where there is a many-body crossover from weak (param- 
agnetic) to strong coupling (screened) as the temperature 
is lowered below the Kondo temperature. In our case, the 
screened state occurs when the impurity spin fully aligns 
with the surrounding background below an energy Ta- 
We are particularly interested in the anomalous relax- 
ation of the Cu ion spins that can be measured by NMR. 
In particular, we calculate the frequency and tempera- 
ture dependence of the imaginary part of the transverse 
impurity susceptibility for the QHAF in d = 3, Xi(w, T). 

At long wavelengths and low energies the magnetic 
problem can be described by spin coherent state path in- 
tegral in terms of the Euclidean action Se (we use units 
such that ?i = fcs = 1) with Se = Swz +Sm where Swz 
is the Wess-Zumino term that describes the quantum dy- 
namics for the impurity spin S and 
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where Xi^ 



T, Xi) with /i = 0, 1, .., c? is the space- 
time coordinate, c = 2\/dJaS is the spin-wave veloc- 
ity, g — (? I Ps is the coupling constant, — JS'^a?^'^ 
is the spin stiffness for the non-linear sigma model [11] 
described by the vector field n (a is the lattice spac- 
ing) , A cx J' is the matrix coupling between the impurity 
spin and the spin environment. The action (1) has to 
be supplemented by the local constraint Yi?{x^) = 1. In 
the ordered phase, we can write n « ((^i,(^2,l) with 
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\^Pa\ ^ 1- We will assume A = (Ai,A2,A) and that 
A, Ai,2 <^ Do J where Dq is the bare cut-off of the 
problem. The action reduces to: 



A53 + ^ A 
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a=l,2 



(2) 



where <p = {(pi,ip2) represents the two Goldstone modes 
of the antifcrromagnct. Eq.(2) describes a problem of 
two free bosonic modes coupled to an impurity via its 
different spin components. Notice that the A coupling 
describes the molecular Weiss field applied by the anti- 
ferromagnet. The other two terms have a very different 
meaning, they represent the quantum fluctuations of the 
impurity due to the coupling to the Goldstone modes. 
Since the operators Sa obey the spin algebra they affect 
the impurity spin by inducing transitions between the 
eigenstates of S^. The problem can be reduced to a one- 
dimension (ID) problem using an s-wave expansion [1] 
with a Hamiltonian, H = Hq + Hi5{x), where: 
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Hi = AS3 + ^-^^Xa dkk^"-'^/^ '^>a{k)Sa , (3) 

where ^a{x) and Ila{x) with a = 1,2 are conjugate real 
scalar fields on the half line ($a (fc) is the Fourier trans- 
form of ^a{x) where k is the momentum). It is convenient 
to further decomposed the fields into right, ^a,R{x), and 
left, ^a.hix), moving components associated with out- 
going and incoming waves out of the impurity, respec- 
tively. Notice that the coupling of the impurity spin to 
the bosonic reservoir depends on the dimensionality. This 
should be contrasted with the Kondo problem where the 
coupling to fermions with a Fermi surface make the prob- 
lem insensitive to d. To make contact with the Kondo 
problem we return to the path integral language and 
trace over the bosonic modes in order to obtain the ef- 
fective action Sgf / = Swz + <5/ for the impurity alone, 
where: 



interactions in a magnetic field. There is a Kosterlizt- 
Thouless (KT) phase transition at 71 = 1: for 71 < 1, A 
scales to infinity indicating that the impurity aligns with 
the bulk - in the Kondo language this is equivalent to the 
formation of the Kondo singlet; for 71 > 1, A is irrelevant 
under the RG and scales to zero - this is the equivalent 
of the Kondo problem with ferromagnetic coupling. The 
second case is 1 < a < 2 and 71 = 72 = 7- Notice, 
on the one hand, that in terms of its Fourier transform 
the second term in (4) behaves like IcjI""""^ |S'a(ct;)p where 
(jj is the frequency. On the other hand, Swz describes 
the area of the unit sphere bounded by the trajectory 
parameterized by S(t). For a variation of S by 5S the 
variation in this term is simply SSwz = SS ■ {S x Ot-S) 
and therefore Swz scales like w. Thus, for a < 2 (d < 3) 
the long range interaction is relevant at low energies and 
one can disregard Swz, that is, the impurity behaves 
classically. It may be surprising that as one lower the 
dimensionality (decreases a) the magnetic impurity be- 
haves classically but this results from the fact that the 
interactions in imaginary time become longer ranged. By 
disregarding Swz the problem reduces to a classical XY 
chain with long-range interactions [12] in the presence of 
a field in the Z direction (proportional to A) where 7 
plays the role of the inverse of the temperature. A is 
a relevant perturbation and the spin always orders with 
the bulk without quantum effects. For the d = 2 QHAF 
this problem has been studied via spin-wave T-matrix 
scattering [13], field theoretical methods [10], and quan- 
tum Monte Carlo [14]. Finally, when a > 2 {d > 3) the 
second term in (4) is irrelevant and the spin effectively 
decouples from the fluctuations of the environment at low 
energies, that is, the problem can be described in terms 
of a quantum spin in the Weiss molecular field alone. 

It is clear that the case of most interest is when a = 2 
and 71,72 7^ 0. Notice that for an antiferromagnet this 
implies d = 3 which is also the case of experimental in- 
terest. Returning to (3), it is convenient to work with 
a single field on the entire real line by the unfolding 
transformation: (I)l{x) = ^i^l{x), ipnix) = '^2.Rix), and 
<I>l{-x) = '^i,r{x), (j)R{-x) = ^2,l{x), all for a; > 0. 
In this case the Hamiltonian can be written as (we set 
c=l): 

+ CSO 



/-TOO 
dx {d,Mx)f + S{x)[AS: 

J ov. ^ - J -• \T._T.,\a ' V- -°° c=R,L 



a=l,2 

where 7a = {XlgSdTid- l))/(4(27r)'^c'') {Sd is the area of 
the hypersphere in d dimensions and r(a;) is the Gamma 
function) and a = d — 1 (for a ferromagnet a = d/2). 

There is a few well understood limits of this problem. 
Consider first the case where 71 ^ but 72 = and a = 2 
{d = 3). In this case the action (4) can be mapped onto 
the anisotropic Kondo problem [3] and is equivalent to 
the problem of a classical ID spin chain with long range 



%1TKldx(t)R{x = Q)Sl- vB'KK2dy(j)L{x = 0)82 , (5) 

where 4ir{x) and 4il{x) are right and left moving fields, 
and Ki^2 = ff^^^Ai_2/(47r''/^) are the new couplings. Simi- 
lar problems to the one described by (5) have been stud- 
ied in the past. In one case a single heat bath is coupled 
to different spin components but in the absence of the 
Weiss field [6,15], in another, the Weiss field is consid- 
ered only in first order [16]. In our case it is not only im- 
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portant to have two independent heat baths but we also 
consider the strong renormahzation of the Weiss field. 

We have performed RG calculations for (5) in the 
Coulomb gas formulation [17] in two different limits: (i) 
1^1.2 ^ 1; Ki 1 with K2 arbitrary and vice-versa. 
The RG equations are: 

diKl = —Kik\ — K1K3 , 

dtK2 = —H2n\ — K2K3 , 

9^K3 = (1 - - K2)k3 , (6) 

where K3 = A/ZJq < 1 and ^ = \ti{D/Dq). Notice that, 
as expected, these equations are symmetric under the ex- 
change of couplings Ki and K2. There are 3 fixed points: 
(1) the trivial one at ki — K2 — K3 = 0, (2) ki = K3 = 0, 
K2 — I, and (3) K2 = ^3 = 0, ki = 1. The two non-trivial 
fixed points are the KT transitions of the anisotropic 
Kondo model. In the ki x K2 plane (K3 7^ 0) the RG 
flow is shown in Fig.l. 
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FIG. 1. RG flow in the m x K2 plane. 

Notice that when the couplings are large and 
anisotropic {ki ^ K2 > 1) the RG indicates that one 
of the couplings flow a fixed value while the others flow 
to zero (k3 ^0). In terms of (5) this indicates that 
the impurity spin aligns in a direction perpendicular to 
the molecular field which was assumed to point in the Z- 
direction in (2). In the Kondo language this is the equiv- 
alent to the Kondo effect with ferromagnetic coupling, 
when the impurity decouples from the environment. A 
possibility that is not considered in this work is associ- 
ated with the formation of a spin texture around the im- 
purity spin. In a classical spin system a spin texture can 
be formed in the bulk spins due to the presence of strong 
and/or anisotropic interactions. The spin texture can 
follow the impurity as it tunnels invalidating the meth- 
ods used here (an instanton calculation is required to take 
into account the collective nature of the texture) [8] . Our 
results are only valid if no spin texture is formed around 
the magnetic impurity. 

In the isotropic case when ki = K2 = n the RG equa- 
tions become: 



n 32 

OlK = —K — KK^ , 

di^3 = (1 - 2k2)^3 . (7) 

Observe that contrary to the KT transition K^^i) always 
scales towards strong coupling indicating the relevance of 
the molecular field (although it decreases initially under 
the RG if K > l/-\/2). However, the RG breaks down 
at a scale i* = ln(Do/^yi) when K3{i*) ~ 1. T4 is the 
crossover energy scale from weak to strong coupling (the 
equivalent of the Kondo temperature). It is easy to see 
that the value of Ta depends strongly on the bare value 
of k{£ = 0). If k(0) < K3(0) the term in (7) does 
not play a role, the flow is essentially the same as the 
usual KT flow and Ta ^ Do[k3(0)]^/(^"^'"'(°^^ « A[l - 
2^2(0) ln(Z)o/A)]. If, on the other hand, k(0) > ^3(0) 
then the k'^ term dominates and K3(£) flows to strong 
coupling leading to: Ta « A(l + 2^2(0) ln(£)o/A))-i. 
We immediately notice that the term in the RG de- 
stroys the KT transition. Unlike the Kondo problem the 
system retains coherence even at large coupling and is 
never overdamped. This is a quantum mechanical effect 
and comes from the fact that the spin operators do not 
commute. While the Si operator in (5) wants to orient 
the impurity spin in its direction, the same happens for 
the 5*2 operator. In a classical system (large S) the spin 
would orient in a finite angle in the XY plane. However 
for a S" = 1/2 impurity this is not possible and the im- 
purity coupling is effectively quantum frustrated reducing 
the effective coupling to the environment. Another inter- 
esting feature of the RG flow is that for 1 > k(0) > l/\/2 
the value of k* = k{1*) « ln(Do/lA)/2 is essentially in- 
dependent of k(0) at energy scale T4. While Ta gives the 
crossover energy scale between weak and strong coupling, 
K* provides information about the dissipation rate, t~^, 
of the impurity dynamics. Our results indicate that for 
k(0) sufficiently large, is independent of the initial 
coupling to the bosonic baths. 
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FIG. 2. x±{'^j ct) ^ function of uo/Ta- 

In order to investigate the dynamical correlations we 
study the frequency dependent impurity spin correlation 
function, xl(w) = {^^ dte"^\[Si[t), Si[0)])] . One 
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of the great advantages of writing the problem in the 
form of an impurity problem as in (5) is that it can be 
accurately studied by numerical renormalization group 
methods (NRG) [18]. In order to perform the NRG cal- 
culation we transform (5) into a ID fermionic problem 
and map the bosonic couplings onto fermionic coupling 
by studying the spectrum of both problems [8]. In Fig. 2 
we plot Xj^(a;)/ijj for different values of 71 = 72 = 7 as 
a function of uj/Ta- We see that the curves collapse for 
7 > 0.4 (quasi-scaling) while deviations are observed for 
small enough 7. This result agrees with the RG since the 
width of x_l{^)/^ is exactly the dissipation rate, r~^, 
which becomes independent of 7 (or k(0)) for 7 > 0.4. 
In Fig. 3 we compare the behavior of the correlation func- 
tion in the isotropic case (71 — ^2 — 0-59) to the behav- 
ior in the strongly anisotropic case, that is, in the DTLS 
(71 — 0.59 and 72 = 0). The differences are striking. 
While in the isotropic case the peak in the response at 
LD — Ta remains, that is the system is underdamped, it 
has disappeared in the anisotropic case where the dy- 
namics is overdamped. 
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FIG. 3. x±i^)/'^ a function oiu/D. 

Since the RG indicates that the transverse couplings 
of the impurity to the environment always flow to k ^ 1 
one could use perturbation theory to calculate Xa.{^)- 
However, perturbation theory in k only generates a Dirac 
delta peak at w = A, that is, = . In order to get a 
finite one needs a non-perturbative calculation. The 
RPA with the bare parameters replaced by the renormal- 
ized ones gives [8]: 

xl(^) ^ (^/2)[arctan(Ar)]-iT^/r 
w (a;2_(j^^)2_ 1/^2)2 _^ 4^2/^2' 

where r^^ « {k*YTa. Notice that (8) reduces to 
a Dirac delta function at w = A as k(0) 0, as 
expected. Wc find that this approximation is good 
for cj <C Ta and also describes well the NRG re- 
sults for all w < i^o when Dq > Ta Dqk*. In 
the zero frequency limit (8) reduces to x±(^ — 0) ^ 
{K*fuj/{TAf + 0{{k*)^) and the Kramers-Kronig rela- 
tion immediately leads to X-l('^ = 0,T = 0) = 7r/[8(l -I- 



(K*)'*)arctan(l/(K;*)2)]1/TA « 1/ {ATa) + 0{[K*f). For 
Dq > uj '> Ta > Dqk* (8) agrees with the NRG results 
giving x_l(^) In the case where Ta ^ Dqk*, in 

the frequency and temperature range: Ta 'ti lu,T 't^ Dq, 
we find: xl(^) ~ TT/[8K^{0)uj\n^{Do/u;)] and x±{T) « 
l/[8K:2(0)rin(i^o/T)] [8]. 

In summary, we have studied a problem of a spin 1/2 
quantum impurity coupled to the Goldstone modes of a 
magnetically ordered system and found that the prob- 
lem maps into a generalization of the DTLS that shows 
no decoherence even in strong coupling. We have cal- 
culated the frequency and temperature behavior of the 
impurity susceptibility that can be measured directly in 
a NMR experiment. We assign the destruction of deco- 
herence (and the KT transition) to a quantum frustration 
between non-commuting spin operators. This result may 
have implications in quantum computation where deco- 
herence effects are detrimental and the use of quantum 
frustration may be explored as a way to avoid decoher- 
ence. 
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